New formulation of 7Z box corrections to the weak charge of the proton 
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We present a new formulation of one of the major radiative corrections to the weak charge of 
the proton - that arising from the axial- vector hadron part of the 7Z box diagram, 5Re 0^ z . This 
formulation, based on dispersion relations, relates the 7Z contributions to moments of the F^ z 
interference structure function. It has a clear connection to the pioneering work of Marciano and 
Sirlin, and enables a systematic approach to improved numerical precision. Using currently available 
data, the total correction from all intermediate states is 5Re 0^ z — 0.0044(4) at zero energy, which 
shifts the theoretical estimate of the proton weak charge from 0.0713(8) to 0.0705(8). The energy 
dependence of this result, which is vital for interpreting the Q WC ak experiment, is also determined. 



As modern parity-violating (PV) experiments press to 
ever improving levels of precision, they remain a vital 
complement to direct tests of the Standard Model at the 
high energy frontier. The classic example of this, involv- 
ing precise measurements of parity violation in atoms, 
led to a remarkably accurate determination of sin 2 9yy . A 
complementary PV electron-proton scattering measure- 
ment underway by the Q W oak Collaboration |l( at Jef- 
ferson Lab has the potential to increase the mass scale 
associated with new physics to 2 TeV or higher, provided 
that the critical radiative corrections are under control. 
In this Letter we present a new formulation of the im- 
portant 7Z radiative corrections which allows for their 
controlled, systematic evaluation. 

Including clcctrowcak radiative corrections, the proton 
weak charge is defined, at zero electron energy E and zero 
momentum transfer, as Q 



Q 



(1 + Ap + A e )(l -4sin 2 6 W (0) + A' e ) 
+ a ww + a zz + a lZ (o), 



(1) 



where sin 2 6 W (0) is the weak mixing angle at zero mo- 
mentum, and the corrections Ap, A e and A' e are given in 
and references therein. The contributions D^f and 
O zz arise from the WW and ZZ box and crossed-box 
diagrams, and can be computed perturbatively. They are 
expected to be energy independent for electron scatter- 
ing in the GeV range. By contrast, the jZ interference 
correction O jZ [E) depends on physics at both short and 
long distance scales. 

In the classic work of Marciano and Sirlin (MS) Q, 
□ 7 z(0) was evaluated in a quark model- inspired loop cal- 
culation using either a "pcrturbative" (p) or a "nonper- 
turbative" (np) ansatz, 



A M z) 2~ B p(np) 



(2) 



where u e (M|) = (1 - 4s 2 ) 
0.23116 in the MS scheme [ 
The perturbative ansatz \ 

B n = In 



and 



Ml 



(Ml) 



is the free quark model result, with m a hadronic mass 
scale, and shows the leading-log behavior. For the non- 
perturbative ansatz, B np = K m + L m is the sum of a 
long-distance part, L r , 



and a short-distance part, K m , 



with 



K m = 



da 



u{l + u/M 2 z ) 



1 - 



a s (u) 



(4) 



Here m is a mass scale representing the onset of asymp- 
totic behavior at large loop momenta, and the factor 
(1 — q s (w)/7t) is the lowest-order correction induced by 
the strong interactions. In Ref. Q L m is taken to be the 
elastic nucleon (Born) contribution, which is evaluated 
to be 2.04 using the same dipole form factors for both 
the electromagnetic and axial-vector coupling. MS Q 
originally adopted the value K m = 9.6 ± 1, based on cal- 
culations with m in the range 0.3-1.0 GeV. A more recent 
calculation by Bardin et al. Q sets 0.5 < m < 0.6 GeV, 
over which K m varies from 9.20 to 9.17 using a 3- loop 
evaluation of a s . Marciano gives an updated value 
for B np of 11.0 ± 1.0, but in view of the high momen- 
tum scales in Eq. (j4}, suggests replacing a by a(M z ) in 
Eq. ^ . This value for U lZ is the one adopted in Ref. [IJ , 
and contributes almost half of the error in the theoretical 
estimate Q p w = 0.0713(8). 

To progress in a systematic way beyond the approach 
of MS and to determine the dependence on energy E, 
we present a new formulation of the box diagram contri- 
bution in which the dominant part of the correction is ex- 
pressed in terms of empirical moments of structure func- 
tions. At forward angles one can compute ^-yz(E) from 
its imaginary part using dispersion relations [7[. The 
imaginary part depends on the PV ep — > eX cross sec- 
tion, which can be expressed in terms of the product of 
leptonic and hadronic tensors. Following standard con- 
ventions [3| , the hadronic tensor can be written in terms 
of the interference electroweak structure functions as 



.fivXp -PAgp pjZ 



p ■ q 



lp-q 



(5) 



(3) 

where p and q are the four-momenta of the proton and 



2 



exchanged boson, respectively. The F^ z contributions 
to d-yz involve the vector hadron coupling of the Z , and 
were recently computed in Refs. [7Hl(j|. 

Our focus here is on the F^ z contribution involving 
the axial- vector hadron coupling of the Z. Following an 
analogous derivation in Ref. H, we can write 



1 



{2MEf 



dQ 1 



v e (Q 2 )a(Q 2 )F. 
Q 2 /M 2 Z 



■yZ 
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W 2 -M 2 + Q 2 
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2 1,(6) 



elastic 



resonance 



DIS (Q <1), Model 1 
DIS (g 2 <l), Model 2 
DIS (2 2 >1, n > 3) 



with s = M 2 + 2ME and Q 2 nax = 2ME{1 - W 2 /s). The 
real part is determined from the dispersion relation 

E' 



5fte d^ z (E) = - dE '-^?2 ^2 Sm U izi E ')i ( 7 ) dashed) contributions. For the DIS part, the high-Q 
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which accounts for both the box and crossed-box terms. 
Unlike the vector hadronic correction die0^ z (E), which 
vanishes at E = 0, the axial- vector hadronic correction 
5fte remains finite, and is dominant in atomic par- 



11 



ity violation at very low electron energies 

We incorporate one further improvement over earlier 
calculations by allowing for the Q 2 dependence of a(Q 2 ) 
and sin 2 9 W (Q 2 ) = k(Q 2 ) s 2 in Eq. ^ due to boson self- 
energy contributions. Both quantities vary significantly 
over the range of Q 2 relevant to these integrals. The pho- 
ton vacuum polarization expression is well-known, and 
expressions for the universal fermion and boson contribu- 
tions to k(Q 2 ) are given in Ref. [U]. Following Ref. Q, 
we use effective quark masses to reproduce the hadronic 
contribution of Aai^Mf) = 0.02786 obtained from dis- 
persion relations [4j, yielding k(0) = 1.030. This is suf- 
ficiently accurate for the purpose of calculating the box 
contributions. In the numerical results that follow, the 
effect of using a(Q 2 ) and v e (Q 2 ) reduces the total contri- 
bution to Eq. ([7]) by 17% relative to using a and v e (M z ). 

The imaginary part of can be split into three re- 
gions: (i) elastic (el) with W 2 = M 2 ; (ii) resonances (res) 
with (M + m 7I ) 2 < W 2 < 4 GcV 2 ; and (iii) deep inelastic 
(DIS), with W 2 > 4 GcV 2 . Contributions from region 
(i) can be written in terms of the elastic form factors as 



F. 



7 Z(cl) 



(Q 2 ) = ~Q 2 G P M (Q 2 )G Z A (Q 2 )S(W 2 - M 2 ). (8) 



For the proton magnetic form factor G P M we use the re- 
cent parametrization from Ref. [ToT ] (the results are sim- 
ilar if one uses a dipole with mass 0.84 GeV), and take 
the axial-vector form factor to be G Z A (Q 2 ) = —1.267/(1 + 
Q 2 /M 2 A ) 2 with M A = 1.0 GeV. A virtue of the dipole 
forms is that the integrals ([6]) and ([7]) can be performed 
analytically, which provides a useful cross-check. 

To simplify notation in what follows, we denote 5fte 
by since that is the quantity of interest in Eq. ([TJ. 
The result for the elastic contribution ( E ) is shown 
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FIG. 1: Real part of D^ Z (E) as a function of incident electron 
energy E. Shown are the elastic (solid) and resonance (dot- 

n > 3 

term (dotted) is neglig ibly small. The two Q 2 < 1 GeV 2 
estimates (long and short dashes) show a very mild E depen- 
dence. Not shown is the dominant high-Q 2 , n = 1 moment, 
which is 32.8 x 10 -4 , and is independent of E. 



in Fig. Q] It agrees exactly with the direct loop calcula- 
tions of Dfz in Refs. 0, Il5j , in which the intermediate 
nucleon is off-shell. It also agrees exactly at E = with 
the value L m = 2.04 if theparameters are adjusted to 
correspond to those of MS [3j. 

For the resonance contributions □^|. res -' from region 
(ii), we use the parametrizations of the transition form 
factors from Lalakulich et al. [l6j], but with modified 
isospin factors appropriate to jZ. These form factors 
have been fitted to the Jefferson Lab pion clcctroproduc- 
tion data (vector part) and pion production data in v 
and v scattering at ANL, BNL and Serpukhov (axial- 
vector part). The parametrizations include the lowest 
four spin-1/2 and 3/2 states in the first and second res- 
onance regions, up to Q 2 = 3.5 GeV 2 . At larger Q 2 the 
resonance contributions are suppressed by the Q 2 depen- 
dence of the transition form factors, which is stronger 
for the dominant A(1232) resonance than for the higher- 
mass resonances [16l |. The resulting resonance contribu- 
tion □^ 1Cb ' > (0) is smaller than the elastic term at E = 0, 
but decreases less rapidly with increasing energy. Vary- 
ing the Q 2 dependence of the poorly determined axial- 
vector form factors has a negligible effect on these results. 

To compute the DIS contributions from region (iii) it 
is convenient to interchange the order of integration in 
(J6j) and (|7|) , in which case the integral over energy can be 
performed analytically A further change of variable 
from W 2 to Bjorken x = Q 2 /(W 2 - M 2 + Q 2 ) gives 



n A(DIS) 



f(r,t) 



2 r 

(E) = - dQ 

k Jo Q 

^dx Fj z (x,Q 2 )f(r,t), 



v e {Q 2 )a{Q 2 ) 
Q 2 (l + Q 2 /M 2 ) 



(9) 



f 2 



[log (1 - i 2 /r 2 ) + 2ttanh~ 1 (t/i 
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with r = 1 + y/l + 4M 2 x 2 /Q 2 , * = 4M#x/Q 2 , and 
^max = Q a /(W£ta - Af 2 + Q 2 ). For t = 0, wc find 



/(r, 0) = (2r — l)/r . In the free quark model limit with 
F 3 7Z = (5/3) x 5(1 - x), Eq. (HJ) then gives exactly the 
perturbative result of Eq. Q for E = (ignoring the Q 2 
dependence of a and v e ). 

To proceed, we divide the Q 2 integral of the full ex- 
pression (j9j) into a low-Q 2 part, where the structure 
function F^ is relatively unknown, and a high-Q 2 part 
(Q 2 > Q 2 ), where at leading order (LO) the structure 
functions can be expressed in terms of valence quark dis- 
tributions q v = q — q Q , 



E 



7 Z(DIS) 



(x, 



%e q g q A q v (x,Q 2 



(10) 



At high Q 2 and low E, the integrand in © can be ex- 
panded in powers of x 2 /Q 2 , yielding a series whose coeffi- 
cients are structure function moments of increasing rank, 



3 
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dQ 2 



M^(Q 2 



2M" 



(5E 2 -3Q 2 )M^>(Q 



v e (Q 2 )a(Q 2 ) 
r(l + Q 2 /M 2 ) 

(3), 



.(11a) 



For completeness, we also quote the result for the vector 
hadronic correction, 



V(DIS) 



(E) = 



2ME 



dQ 2 



a(Q 2 ) 



Q 4 (l + Q 2 /M 2 ) 



M< 2 > (Q 2 ) + \m[ 2 Xq*) + ^(E 2 ~ Q 2 )M^\Q 2 ) 



2M 2 
~~5Q' 



-(4E 2 -5Q 2 )M 1 (4) (Q 2 ) + 



(lib) 



In Eqs. ([lip the moments of the structure functions are 
defined as 



Ml n \Q 2 )= / dxx 
Jo 



i-2 -rnZ 



1,2,3, (12) 



where Tj Z = ^xF^ z , F^ Z ,xF^ z ^. In approximating 

the upper limit x max on the x-intcgrals in Eqs. (|11[) by 
1, the resulting error is less than 10~ 4 for Q 2 > 1 GcV 2 . 
The large- x contributions to M^ n \Q 2 ) become more im- 
portant for large n; however, the higher moments are 
suppressed by increasing powers of 1/ Q 2 . In practice, the 
integrals in Eqs. (TTTj) are dominated by the lowest mo- 
ments, with the 1/Q 2 corrections being relatively small 
in DIS kinematics. 

Equations (jlip are major new results which provide a 
systematic framework within which to evaluate the ra- 
diative corrections. For the axial- vector hadron part, the 
lowest moment, M^\q 2 ), is the analog of the GLS 



sum rule [17[ for vN DIS, which at LO counts the num- 
ber of valence quarks in the nucleon. The corresponding 
quantity for 7Z is J2 q ^ e q 9 a = so that at next-to- 
lcading order (NLO) in the MS scheme 



M«(Q 2 ) = \(l a ° {Q2) 



(13) 



Mf (O 2 



)= i( 2( ^ +( x 2 ) d )(l + ^)), 

where (x 2 ) q = J Q dxx 2 q v (x, Q 2 )- Hence the lowest (n = 
1) moment contribution to Eq. (|lla|) is identical to the 
MS result Q in Eq. ((4]). However, the parameter Qq in 
Eq. (jlla.p has a slightly different interpretation than the 
mass parameter m 2 of Eq. (U). Here Qo corresponds to 
the momentum above which a partonic representation of 
the non-resonant structure functions is valid, and above 
which the Q 2 evolution of parton distribution functions 
(PDFs) via the Q 2 evolution equations is applicable. We 
take Q 2 , = 1 GeV 2 , which coincides with the typical lower 
limit of recent sets of PDFs 0, [l!| . The computation of 
the vector hadronic contribution to □^ IS ' ) proceeds in a 
similar manner, and will be discussed elsewhere poj ]. 

To evaluate the moments in Eq. (jllap we use several 
NLO parametrizations of PDFs determined from global 



fits [18|,[19j]. The results are summarized in Fig.[TJ Varia- 
tions in the values of a s (Af§) among the datasets consid- 
ered had a negligible effect on the n — 1 value of 0.0033. 
The n = 3 moments for different datasets are virtually 
identical, and give negligibly small contributions. 

The E dependent terms in Eq. (jllap should also be 
small, since these depend on n > 3 moments. However, 
the expansion in Eq. (jllap is not strictly valid when E > 
Qq/2M. To describe the E dependence in this region we 
evaluate the difference □^ DIS) ( J B)-D^ DIS) (0) in Eq. © 
by replacing f(r,t) by f{r,t) — f(r,0). The results are 
indeed small for E in the few GeV region, as the dotted 
line in Fig. [T] indicates. 

For Q 2 < Qq a partonic description of the structure 
functions is not valid. In particular, since the integral 
over Q 2 in Eq. (j9)) extends down to Q 2 = 0, and the upper 
limit on the x-integral, x max , is also limited by Q 2 , one 
requires the behavior of the structure functions at both 
low x and low Q 2 . In the case of the vector F^ Z structure 
function, conservation of the two vector currents requires 



F^ 2 ~ Q 2 as Q 2 — > 0. By contrast, F^ depends on both 
vector and axial-vector currents, and the nonconservation 
of the latter means that no similar constraint exists [lij ]. 

In the absence of data on F^' z (x,Q 2 ) in the low-x, 
low-Q 2 region, we consider models for the possible x and 
Q 2 dependence, obeying the following conditions: (1) 
,Q 2 ) should not diverge in the limit Q 2 — > 0; 



F lZ (r 

3 V^niax? 



(2) F^ z (x, Q 2 ) should match the partonic structure func- 



tion at Q = Qq. For the paramctrization of Ref. [Uf we 



note that F^ z (x, Q 2 



-0.7 



as x — > 0. With this in 
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FIG. 2: Total (el+res+DIS) axial-vector hadron correction 
a^ z (E) (labeled "A") and the sum of axial and vector hadron 
[8] corrections (labeled "V+A"), together with the E = 
result of MS || (extended to finite E for comparison). The 
vertical dashed line indicates the energy at Q wca k kinematics. 

mind, we consider two models for Q 2 < Qq. 



Model 1 sets 



Ff{x,Q 2 



1 + A 2 /Q 2 



(14) 



which has the property that F^ z (x n 



l ) - (Q ) 



2\0.3 



as 

Q 2 — > 0. Here A 2 is a parameter that can be adjusted to 
examine the model sensitivity of the integral in Eq. © . 
For A 2 in the range (0.4 - 1.0) GeV 2 , we obtain a ±10% 
variation in the values for n^ z (E) shown in Fig.Q] 

Model 2 freezes F^ z at the Q 2 = Q\ value for all W 2 , 



F^(x ,Q 2 ), 



which is equivalent to setting F^ z (x, Q 2 ) 
with x = xQl/ ((1 - x)Q 2 + xQl). For this model, K 
is constant as Q 2 — > 0, and yields a 15% larger contribu- 
tion to U^ Z (E) than Model 1, as illustrated in Fig. [U 

The total correction to is given by the sum 

(el+res+DIS), and is shown in Fig. Has a function of 
E. As demonstrated, the E dependence arises predomi- 
nantly from the elastic and resonance contributions. We 
assign a very conservative uncertainty estimate equal to 
twice the low-Q 2 DIS value. This allows for uncertainties 
in the resonance and low-Q 2 DIS contributions, and in 
the effect of the running coupling constants on the dom- 
inant n = 1 contribution. The total contribution to t3~ z 
is 0.0044(4) at E = 0, and 0.0037(4) at E = 1.165 GcV 
(the Qweak energy). This should be compared to the 
value 0.0052(5) used in Ref. 0, which is assumed to be 
energy independent. Also shown in Fig. [2] is the total 



J 7 Z 



using the result for from Ref. 



which has an uncertainty that grows with E. 

Our value shifts the theoretical estimate for Q 1 ^ from 
0.0713(8) to 0.0705(8), with a total energy dependent 
correction U lZ (E) - Oyz(0) of 0.0040]^;^ at E = 
1.165 GeV. A similar uncertainty would be obtained us- 
ing the estimate of tiY z from Ref. Q, while a larger un- 
certainty on the vector hadron correction was quoted in 



Ref. [10( . These uncertainties can be reduced with future 
PV structure function measurements at low Q 2 , such as 
those planned at Jefferson Lab. The high precision deter- 
mination of would then allow more robust extraction 
of signals for new physics beyond the Standard Model. 
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